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ABSTRACT 

The finite element method is applied to the vibration 
analysis of isotropic and polar orthotropic plates. Annular 
finite elements having four degrees of freedom, these being 
the displacements and slopes at inner and outer radii, are 
used to obtain the natural frequencies of all the nine combi- 
nations of clamped, free and simply supported edge conditions 
of annular plates. Using the annular elements itself, the 
natural frequencies are obtained for solid circular plates 
having clamped, free or simply supported boundaries. 



CHAPTER 1 . 


INTRODUCTION 

1.1 GENERAL 

Vibration characteristics of circular plates are of 
interest in design since circular plates with and without 
holes are commonly used structural elements. Fibre-reinf orced 
composite materials are attractive for these structural ele- 
ments because of their high stiffness, relatively low density 
and highly controllable anisotropic properties. Moreover 
polar orthotropic composites have great potential for the 
development of flywheels for various applications. The 
transverse flexural vibration can constitute a major problem 
in designing flywheels. Analysis of these polar orthotropic 
materials generally defies the possibility of exact solution 
and necessitates the use of approximate methods. Exact solu- 
tions are possible only for highly restricted classes of 
problems? the restrictions tend to be more severe for vibra- 
tion problems than for static problems. 

In the last 20 years high speed electronic digital compu- 
ters and finite element methods suitable for computer aided 
analysis have been rapidly developed. The governing differen- 
tial equations describing the behaviour of any continuum are 
derived considering infinitesimal elements of the continuum. 



2 


A finite number of these infinitesimal elements build up the 
continuum which thus has infinite degrees of freedom. In 
finite element idealisation, elements of finite size are 
considered in place of infinitesimal elements to build up 
the actual continuum, and the behaviour of the continuum is 
sought to be represented through the behaviour of individual 
elements. In this idealisation each element is connected to 
its neighbouring elements, not continuously, but through a 
finite number of strategically located points or nodes on its 
boundary, thereby reducing the infinite degreesof freedom of 
the continuum to a finite number of degrees of freedom. Using 
the same principles that also hold good for the actual conti- 
nuum, the behaviour of each element is exactly or approximately 
determined in terms of the variables at the connecting nodes. 

If the required continuity of the variables is maintained bet- 
ween the elements, then the finite element solution is expected 
to converge towards the true solution with greater number of 
smaller elements. So, in order to study the dynamic response 
characteristics of polar orthotropic plates the finite element 
method is found to be less tedious compared to other approximate 
analyses and hence in this present work finite element method 
is used for the analysis. 

1.2 REVIEW OF PREVIOUS WORK 

In the literature, one finds an extensive treatment of the 
problem of fluxural vibration of isotropic circular plates and. 



to a lesser extent, annular plate as well. The problem of 
vibration of isotropic circular plates having free edges was 
investigated long back by Kirchhoff, Lamb and Rayleigh as cited 
[l ], using Poisson-Kirchhoff theory. Timoshenko [2 3 used the 
energy method for solving the case of isotropic plate with 
clamped edges, Thein Wah [l] investigated the case of simply 
supported edge taking Poisson-Kirchhoff theory as the basic 
equation. Numerical results are also given for frequency 
parameter for modes of vibration consisting of nodal circles 
and nodal diameters. An extensive study of uniform annular 
isotropic plates has been done by Vogel et al [s]. Here the 
classical theory of fluxural motions of elastic plates was used 
to determine the natural frequencies for all nine possible com- 
binations of clamped, free and simply supported edge conditions 
and for different hole sizes. An investigation was also made 
into the effects of Poisson's ratio on frequency parameter by 
taking different values of Poisson’s ratio. No significant 
variation was noted in the frequency parameter for different 
values of Poisson’s ratio. Hence, in reference [3 ]» the analy- 
sis has been done by taking Poisson’s ratio equal to 0.3. 

In the text by Zienkiewicz [4 ] several applications of the 
finite element method {FEM) to the vibration analysis of circular 
plates are given. The sinplest element shapes for plate pro- 
blems are obviously a triangle and rectangle with three and four 
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nodes» respectively [5 ]. Rectangular elements are somewhat 
limited in their boundary shape applicability, whereas trian- 
gular elements may be used to represent any boundary. However, 
in our present problem, the use of triangular elements means 
that the curved boundary is being approximated by a series of 
straight line segments. In order to overcome this approximation 
Ergatoudis [ 6 ] has developed quadrilateral elements with 
various curved sides. Here , the basic shape of the element 
chosen is quadrilateral, the sides of which can, however, be 
distorted in a prescribed way, M.D, Olson et al [7 ] have 
developed two plate bending finite elements in polar coordinates 
mainly to study the vibration problem of circular plates. The 
first element has three nodal corners and is in the shape of a 
sector of a circle. The second element has four nodal corners 
and is in the shape of a segment of an annulus. The transverse 
displafcement and two slopes are used as the degrees of freedom 
at each element nodal corner, resulting in nine degrees of 
freedom for circular sector element and twelve degrees of 
freedom for annular sector element. Non-dimensional frequency 
parameter for corqslete circular plate with clamped, free and 
sinply supported boundaries are tabulated and compared with 
exact values. In a recent paper by Kirkhope et al [8 ] 
annular finite elements are derived which describe the bending 
and stretching of thin rotating discs. These elements incorpora- 
te the desired number of diametral nodes in their dynamic 
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deflection function, and allow for any specified thickness 
variation in the radial direction* 

Even though the literature contain a good number of 
treatment for the vibration of isotropic plates, the vibra- 
tion analysis of orthotropic circular plate have, however, 
apparently received scant mention in the literature. Kirmser 
et al. [9 3f Pandalai et al. [lo],Joung [ll], Minkarah ct al 
[ 12 ], Bellini et al. [l3], Huang [14] and Nowinski et al.[l5] 
have attempted some of the problems of polar orthotropic 
bodies subjected to dynamic loading. It is to be noted that 
the existence of material singularity at the centre has been 
reported by many of them. This means that the circumferential 
and radial moduli at the centre cannot be different for polar 
orthotropic materials. In order to include the aforesaid 
material singularity at the centre some investigators have 
assumed an isotropic core of uniform property vdiile dealing 
with the problem of circular disc. 

Natural frequencies corresponding to axisymmetric mode 
of vibration of polar orthotropic annular plates have been 
obtained by Vijayakuraar et al. [16] for various combinations 
of clamped, free and simply supported edge conditions. Classi- 
cal Rayleigh-Ritz method with sin?>le polynomials as admissible 
functions has been used to get the least eigenvalues. In an 
attenpt to determine the natural frequencies of higher modes 
of vibration, Ramaiah et al. [l7] have used the classical 
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Reyleigh-Ritz method by introducing a coordinate transfor- 
mation* Thus they have obtained frequencies corresponding 
to the axisymmotric modes and modes with one and two nodal 
diameters. However, they themselves have admitted that the 
application of this method is generally cumbersome since the 
method involves solutions of many simultaneous equations, 
if reliable estimates of the frequencies are desired. It 
may be due to this difficulty that the numerical results are 
given only for simply supported-free, clamped-free and free- 
free annular plates and this also for two selected ratio of 
radii 0,1 and 0»5, They have determined these frequencies 
for two different cases of rigidity ratio 0.5 and 2.0, In 
another recent investigation by Ramaiah et al, [l8 ] the 
frequencies corresponding to higher modes of vibrating polar 
orthotropic annular plates have been obtained based on an 
assumption that the radial bonding moment is zero at a nodal 
circle. No numerical results are given. 

1.3 PRESENT WORK 

What is lacking in the literature, is a thorough numeri- 
cal evaluation of natural frequencies of polar orthotropic 
circular plates with simply supported, clamped and free edge 
conditions, and frequencies for polar orthotropic annular 
plates with all the nine combinations of simply supported, 
clanged and free edge conditions. To evaluate those natural 
frequencies, in this present work finite element method is 



7 


used which is found to be more efficient when compared with all 
other methods. 

The study of natural frequencies of polar orthotropic 
circular as well as annular plates is the concern of the 
present work. The FEM is applied to the vibration analysis 
and the natural frequencies are computed by means of a 
finite element program. Kirkhope et al. [s ] have recently 
developed annular finite elements to describe the bending 
and stretching of thin rotating isotropic discs. These ele- 
ments allow for any specified thickness variation in the 
radial direction. In the present work this approach has been 
adopted to determine natural frequencies of non-rotating uni- 
form isotropic plates and its applicability is extended to 
polar orthotropic plates also. In order to formulate the 
eigenvalue problems associated with the vibration analyses, the 
element stiffness and mass matrices are derived for isotropic 
as well as orthotropic plates by appropriate manipulation of the 
chosen displacement function. 

In Chapter 2 the derivation of stiffness matrix for iso- 
tropic plate is given. It is derived by putting the assumed 
displacement function in the expression for strain energy. 

The mass matrix is derived by substituting the assumed displace™ 
ment function in the equation for potonti-al energy of the plate. 
Non-dimensional frequency parameters are obtained for all the 
nine possible combinations of clamped, free and simply supported 
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edge conditions of the annular plates* The results are 
compared with exact values and graphs are plotted to find 
the variation of frequency with hole size. The frequencies 
of solid circular plates with clamped, free and simply 
supported edges are also calculated and compared with exact 
values . 

In Chapter 3, the stiffness matrix for polar orthotropic 
annular plate is derived. The expression for mass matrix of 
polar orthotropic plate will be the same as that of isotropic 
plate because the mass matrix is independent of material pro- 
perties Dj.*D0 arid D^. The frequencies of polar orthotropic 
circular plates for all the three edge conditions are determined. 
For polar orthotropic annular plates the frequencies are obtained 
for all the nine cases of edge conditions and for different hole 
sizes. These frequencies are calculated for a wide range of 
rigidity ratio and are compared with the available values. 

Results and discussion are given in Chapter 4 and in 
Chapter 5, some general conclusions are summarized and possible 
extensions of the present work are discussed. 
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CHAPTER 2 


GENERAL FOf^AULATION OF VIBRATION PROBLEivi OF CIRCULAR 

AND AiANULAR ISOTROPIC PLATES 


2.1 INTRODUCTION 

One of the problems in structural is the vibration 

of circular and annular plates. Exact circular[ 1 ] 

and annular plates [ 3 1 for various edge conditions are readily 

^ methods used to study 

available. One of the several approximate 

the vibration problem is finite element method (FEfA) , Olson 
et al [ 7 3 used circular sector elements having nine degrees 
of freedom and annular sector elements having twelve degrees of 


freedom to estimate the natural frequencies 


of isotropic plates. 


These eleirents result in a very large number of degrees of 
freedom if good accuracy is desired. The annular element deri 

ved by Kirkhope et al [s 3 for rotating dis^ of nonuniform pro 

^ j ^ ^ , T+ has been observed 

file has only four degrees of freedom. 1^ “ 

• 4.U * 1 1 4. results obtained are 

that even with two annular elements,- the re& 

1 4. 4.U i I TT the present chapter, 

very close to the exact values. Hence, m ^ 

. , - . - . ^ _ nnd annular plates 

natural frequencies of isotropic circular 

4 . 4 . 4 , - I 4 -- -4 , ^r,ts having four degrees 

are determined using annular finite elemen^i» 


of freedom 
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2.2 FOK^ULATION OF THE ANNULAR ELENENT 

The vibration of circular plates is characterised by 
modes having integer numbers of nodal diameters and nodal 
circles. For each nodal diameter configuration there exists 
a family of circular nodes with corresponding natural fre- 
quencies. Fig. 2,1 shows the annular finite element with 
associated degrees of freedom and diametral nodes. The 
finite element has two nodes, with 2 degrees of freedom/node. 
They are : the axial deflection ( w) , and radial rotation (0). 

The rotation 0 is defined as 0 = and the bars over w and 0 
indicate that these are taken along the reference antinode. 

Once, w, and 0, at the reference antinode, where 0 is taken as 
zero, are specified, the deflection and slope at any other 
point at an angle 0 from the reference antinode, as shown in 
Fig. 2.1, can be expressed as w cosm© and ^ cosm0. Hence, a 
suitable deflection function for w of the plate along the radial 
direction only remains to be chosen. 

The exact solution to the problem is of the form 
w(r,0) = F(r) cos m©. So we can choose a deflection function 
of the element, as w(r,0) = F(r) cos rn0, where F(r) is the 
deflection of the vibrating plate along the reference antinode 
and can be approximated by a polynomial of r. Use of cos 
in the expression incorporates the desired number of nodal 
diameters m in the element. 
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2.2*1 Stiffness Matrix for Isotropic Plates 

For classical bending of isotropic plates, the strain 
energy, for an element in polar co-ordinates is given by [2 ], 

TT „ D r r/d^w x2./l dw . 1 d^w\2 . d^W/ 1 dw . 1 d^w\ 


)(| IjjQ - §g)^]rdrd8 (2.1) 
r 


which can be written as 
2ti a 

U=| / / {dj [V] fdj rdrdO 

0 b 

19 2 2 

/d W\/1 dw ^ 1 d W\/l d w 

where, idj = . ^ ^ 


1 ^w^ 


( 2 . 2 ) 


(2.3) 


0 2{l-i ) 


(2.4) 


where is Poisson’s ratio. 


For the annular element the assumed deflection function 
can be written as, 

w(r,0) = F{r) cos n®, where F(r) is a polynomial of r, 

2 3 

w(r,0) « ( aj|_+a 2 i‘+a 3 r +a^T ) cos mO 

w(r,9) « UJ |c} cos rfl9 (2.5) 
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where , 

LsJ » U- r 

and 

1 ®2 ®3 

Since 0=0 along the reference antinode, the lateral 
deflection and radial slope along this antinode are, 

w(r) = aj^+a2r+a2r +a4r 

^(r) * a2+2a2r+3a4r^ 

Referring to Fig. 2,1 

— 2 3 

* ^'*■^^3^ t3a^b^ 

— 2 3 

^2 * 3^^+323+833 + 3^3 

^2 * 82+2333+33^3^ 

Let ^w| be 

■^1 
^1 
% 

w| » [E] { c j 




W * [a 


then 


( 2 * 6 ) 



14 


where 



From eqn. (2.6) 

(c] = [E]-^ {«] 

(c] » [B] (2.7) 

where [B] =* [E]"*^ 

Putting this value of ic v in eqn. (2.5) 

w(r,0) » ^sj [b] ^^w^cos n© (2.8) 

[b3 = 
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Putting this value of w(r,0) in equation (2,3) 




0 

7" 


0 


2a 


6ra 


~(m^-l) -a(m^-2) -ar(r/-3) 


0 


-mp 


-2mpr 


[B][ 


w 


'j 


where a « cos m9 


p * sin 


(d) = [F] [B]{w} 

Substituting this value of ;^dj in eqn. (2.2) 


(2.9) 


271 a 


u = i ; / LwJ [ef [f]T [V] [F] [B] £w] rdrd9 


o b 


27t a 


- 5 H [B]^ C / / [F]^ [V] [F] rdrdG] [B] [«} 
o b 

Stiffness matrix is given by 

2tt a 

[K] = [b]T [ / /D[Ff [V] [F] rdrdS] [B] 

o b 

[K] » [Bf [k] [B] 


271 a 

where [k] -D f / [f]^ [v] [f] rdrd© . 

o b 


16 . 


2% a 

Ek] ^T>S S 

o b 


1 

i 

0 


2ic a 

=T>/ / 

0 b 


-•2a 


-6ar 


2 

m a 
r 

r 


a(ni^ -2) 


r 

0 

-mp 


ar(ni -3) '-2mpr 


O 

0 


o 


0 

0 

-2a 

-6ar 

m^a 

I (m2-l) 

a(m^-2) 

2 

ar(m -3) 

? 

0 

-mp 

-2inpr 


*^am^ 

■ 


x(rdrcl9 ) 


2 

m a 

7“ 


2(1^ ) H| 


r 


I (n?~l) 0 

-2a+'i a(m^-2) -2*^ a+a{m^--2) -2inp(l-^) 

-6ar4%cr(m^-3) ' -6'tar+ar(m^-3) -4n^r(l-^) 

0 -2a -6ar 


2 

r 


■r' 




0 


a(ra^ -2) ar(m^ -3) 


, -mp 


-2inpr 


rdrdO 
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11 

^12 

^13 

^14 

21 

^22 

^23 

^24 

31 

^32 

^33 

^34 

A1 

^42 

^43 

^44 


2ti a -4^2 ^ 

^11 ^ + 2(1-'1) S^)rdrd0 => 2it if 


zit 

=D/ (niV+2(l-n)mV)[— i; d9 

o -2 ^ 

*I> 1 i,(nv4+2(l-i )m^) 

(ab) 


[k, , - CkT)( n,W-2i m2) 

1 2{abr 

3 ^ 2 . Q rt 

ssl)/ / ~y i!r(m"^**l) rdrd© 

0 b 

kj^2 =3)Cit(m'^-m2) is^ 


for 



2 it a ^^2 2 2 , 4 « 2 v 

1^3 «]>/ / - aii- 

O b i A A 


^ 2.2 
Im 


=» CJcD(~2^£n^+(!a^--2iri^)-2(l--^ )m^) ln(r2/rj^) 

^ ^ ^ . 

^13 ^(®'^'"4m^) ln(a/b) 


m = o 

mm > 1 


2 if m =* 0 
1 if m > 1 


rdrd© 



rdrdG 


In a 2 2 2, 2 . 2„2,, , 

k . =D/ / (-6 + HL2„(in -3) _ 4m~g„( 1- >) ) 

rv K ^ ^ 2? 


=DCn( -•6'^ m^+ ( in^~3ra^ ) ~4m^ +4^*^ ) ( a -b ) 

* CnD(m'^-7m^-2*ilin^)(a-b) | 

2 % a ^2/ 4 2v 
==•])/ 7 rdrde . 

o b T 

»3Cn(m'^-ii?) (fg^) 


*^21 * *^12 


^22 ^ (m^-1)^ rdrd© 

.0 b 


^22 ~^C7t(m'^-2m^4-l) ln(a/b) j 


^23 (m^-l)+ (m^-2) (m^-l)) rdrd9 

0 b 

» Ciili(-2*^ (ra^~l) + (m^~2) (m^-l)) (a-b) 


^23 OffD(in^-3ra^~ 2'^ra^+2"1+2)(a-b) 


krt^ =j) / / (--6a^'^ (m^-l)+a^(m^-l) (m^-3)) rdrd© 

0 b 

2 2 

g OiD(>-6'1(m^-l)-i-Cm^-l)(ia^>3)) (^ -^) 

1^24 ^•6'^+3) (— ■^-~) 
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^ 27C a ^ ^2 2 2 . ^ „ 2 o2 V 

*"31 + 2^. (rn^~2in2)- (l-^n))rdrd9 

»DCu(m'^~4m^) in(a/b) 
j*"31 ’*"*^ 13 " 

2lii a 2 2 

5^32 J (-2 -~ (ni^-l)+ ^ (m^-2)(m'^-l))rdrd0 

o b 

=3Cit(m^-3m^-2 m^+2 +2)(a-b) 


= k^a 

^32 ~ ^23 

k-- =3)/^ ^ (4a^-'2i:a^(m^-2)-2"^a^(ii?-2)+a^(m^-2)^ 
o b 

+ 2mV(i“'^)) rdrd0 

2 2 

=DC'it(4-2"^ (n?-2)-2'^ (n?-2)+(m^-2)^+(l-'^ )2m^) (^-^— ) 

2 2 

k^a =1 >Cjj ( m^~6 'i "2ra^ +3*^ +8 ) (■-'■ ^ - - ■) ; 


k 


34 


*3)/” (12a^r“6'^a^r(ni^-2)-2a^"^r(m^-3)+a\(m^-3)(m^-2) 

o b 

+ 4(l-''*i)inPp^) rdrd9 

3 3 

:*J)OJt(12-6'i (m^-2)-2'i (m^-3)+(m^-3) (m^-2)+4( l-"i )m*^) (^-^) 


3 3 

k^. *DOjt(ni'^--m^-12"l n?+18’^+18)(^-^) 


^34 
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2 % a ^-4 2 2 2 2 ^ . 2„2 ^ 

1=41 =»/ / (- + a-^ (m 2 - 3 )- ia-i- (l-i))rdrde 

o b ■‘' ^ 

«3)CTt ( m'^-7ra^-2'^ ) ( a-b ) 


^41 ^14 

^42 =1^/ J ('•6'^a^(in^-l)+a^(in2-i)(m^-3))rdrde 
o b 

2 2 

-5 Cn ( m‘^-4m^~6'ii +3 ) 

“^24 

[^42 = ^24 I 

2tc ft 

k43 *])/ / (i2a^r-'2^a^r(m^-3)-6^a^r(m^~2)+a^r(m^_--2)(m^~3) 
o b • 

+ 4m^p^r(l-'*^))rdrd0 
3 3 

*:3Cre ( in'^-rn^-12'^ ii?+18"*+18 ) 

* ^34 

^43 ~ ^34 I 

2ic .3 

^44 *1^/ / (36a^r^-6a^r^'^(ni^-3)-6*^a^r^(in^-3)+a^r^(ra^-3)^ 

0 b 

+ ( l-."i ) ) rdrde 

»DCjt(36-6"iia^+18 '^-6‘im^+i8*^+m'^‘-6n?+9+8m— (^— ) 

1^44 ~l>C>(m'^+2m^>20"^m^-l-36"^-f-45)(^— ) j 
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2.2,2 Mass Matrix for Isotropic Plates 

The mass matrix is derived from the given displaceraent 

function. For classical bending of circular plates the p-etnia- 
Kvm-ulT'c. 

■fetefd energy for an element in polar coordinates is given by [2 j, 

ft K ^ • 0 

T / / w r dr d9 

Ob 

where $ is the mass density. 

Equation (2,8) gives the ddisplaceraent function as 

cos m© 
w = LsJ [B] I cos 

T / / (IwJ[b]^ LsJ[b 3 cos^m9)r drd© 

X [B]^ [ / / |sj cos^m© r drd©][B] (w| 

2 .^ 3 

Mass matrix [M] =^ h[B]^ [ / f L^J cos^ra© r drd©][B3 

o b ^ J 

[M] = [B]f [m] [B] 

where 

2-^ a ^ y . . o 

[nv] - yf h / / [s j JsJ cos m© rdrd© 


w(r,9) =s LsJ [b] ^ 
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2^2,3 Determination of Natural Frequencies 

It should be noted that the element stiffness matrix [k] 
and element mass matrix [m] are independent of global boundary 
conditions. Conversions of [k] and [m] into the total stiffness 
matrix [k] and total mass matrix [m] of the plate is a simple 
mechanical operation of addition. All the relevant boundary 
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conditions are introduced by suppressing the corresponding 
rows and columns of [K] and [m]. Then the governing equation 
of motion for undamped free vibrations of multidegree freedom 
system can be derived as 

Cm] + M {w| = 0 (2.10) 

For a periodic motion the time coordinate is separated by 
assuming the displacement vector as 

{w] = {«o) 

where is the amplitude vector ^ u is the natural frequency 

and ’t’ is the time coordinate* Using equation (2.11) in 
equation (2»10), we can obtain the linear homogeneous equation 
of motion as 

(CK]-[M]aj^) = 0 (2.12) 

This is a linear eigenvalue problem as none of the two 
matrices [K] and [m] involve w. Both [K] and [m] are real 
symmetric and positive definite matrices implying real and 
positive eigenvalues, which in our case are the natural 
frequencies* 

Non-dimensional frequency parameters are calculated 
using a conputer programme in FORTRAN language. IMSL sub- 
routine EICSF is used to get all the natural frequencies 
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2 

of the homogeneous equation of motion (K-fAw ) = 0. These 
frequencies are then converted into nondimensional form. 

For this, in the case of isotropic plate the frequency is 
multiplied by (^h/D)^*^ and for orthotropic plate it is 
multiplied by a^( jPh/op*^*^ . 

The results and discussion are given in Chapter 4. 
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CHAPTER 3 

GENERAL FOFMULATIOM OF VIBRATION PROBLEM OF CIRCULAR 
AND ANNULAR POLAR ORTHOTROPIC PLATES 

3.1 INTRODUCTION 

There are quite a few investigations, reported in the 
literature, on vibration of polar orthotropic circular [9-14] 
and annular platesi 16-18] .Most of these investigations are, 
however, limited mainly to axisyiranetric modes, for particular 
set of edge conditions, and in some of them, for some specific 
values of material constants. In one of these investigations 
Ramaiah et al [l8 ] have obtained natural frequencies of higher 
inodes of vibration based on an assumption that the radial 
bending moment is zero at a nodal circle. But in this analysis, 
the presence of a free edge yield inaccurate estimates of 
frequencies of modes with one or two nodal circles. The exis- 
tence of singularity at the centre of a polar orthotropic 
circular disc has been treated in different ways by different 
investigators. Carrier et al [l9 ] have considered an isotropic 
homogeneous core extending upto 0.1 times the outside radius 
at the centre while Sinha et al [20 ] have taken the size of 
core to be 0.01 times outside radius. In the present analysis 
the question of material singularity does not arise because 
here the frequency is determined by considering the disc as 
an annular plate having radius ratio 0.001, Here, a similar 
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annular element formulated in Section 2,2 has been used to 
determine the frequencies of orthotropic plates. 

3.2 STIFFNESS MATRIX FOR ORTHOTROPIC PLATES 


For classical bending of orthotropic plates, the strain 
energy, for an element in polar coordinate is given by [ ] 


U, 


^ ^ Trk / 2 . n /I dw , 1 ^ W s 2 . r> d w /I dw 


|-|) + 2D^(| - ■ij ^)^]rdrd9 

r Or r 


(3.1) 


which can be written as 


, 2n a ^ - 

Uj. = ^ / / kl [Vil |dj 

where L<lJ = ^ 


rdrd© 



and 

[V,] 





O 


0 


2D 


r© 


(3.2) 

(3.3) 

(3.4) 


2 3 

Assuming a deflection function w(r,©) =( aT|^+a 2 r+a 3 r +a 4 r )cosm9 
and preceding with the analysis given in Section 2.2.1, the 
element stiffness matrix for orthotropic plate can be obtained 

2-h: a 

[k] -IS [F]^ [V, 3 [f] rdrd© 
o b 


as 
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[k] 


= TJ 


o b 


0 


m g M 


0 ~ (m^-l) 0 


-2a “(m -2) 

-6ar ar(ra^-3) -2mprl 


r°r 




0 


"9 


0 


0 




mB 

7 


M = 


2it a 

; f 

o b 


-SgD^m « 




-2a 


*©^r r 


-D^2a+lQD^a(nr-2) 
-6Dj.ar+ "^QD^^arCm^-S) 


-mP 


-6ar 


2L« ^(m^-1) a(m2-2) ar(m^-3) 

_2 


-2ir^r 


rdrd© 


2 

m a 


^© 


2Dj^mp 


Dg f{m2-l) 


-2i.D a+O^a ( ii?-2 ) 

2 

-B'^gDrar+DparCm -3) 


0 


-2D^mp 


-4D:j^niprj 


0 

0 

-2a 

-6ar 


n,2„ 

m a 

|(m2-l) 

a(ni^— 2) 

ar(in^-3) 

rdrd© 

4 

0 

-rap 

-2mPr 
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^11 


^12 

^13 

^14 i 

i 


^21 


^22 

^23 

k ! 
^24 

[k] = 









^32 

^33 

^34 


*^41 


^42 

^43 

! 

If ’ 

44 ; 






— J 

where. 

2% 

a 


4 2 

2D o 

^11 ^ 

f 

r\ 

f 

K 

^ D© 

m a 
.4 

r© 

^4 




rdrd0 


= Cti (D^m + SD^qUI )( 


2x/ a^-b^ 


tQ“‘ / ^ ^ o ) 

2(ab)^ 


11 


Cit(m%+2A^)(4:=^) i 

^ ^ 2(ab)^ : 


^12 


2ic a 2 

/ / [Dq (m^-1)] rdrd9 


o b 


kj^2 = C7i(Dg(ni"'-m^))(a^) 


^13 


2it a 

/ / (• 

o b 


2 QD^m^a^ DQa^(in^-2m^) 

~ Y + -"Y" 

r r 

2 2 
2D 

- — — -j- — )rdrd0 


I kj ^3 = Cre(-2 gD^m^+Dg(m'^-2m^)-2D^ti?) Xn(a/b) j 
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^14 * 


a r-. 2 2 ^22/2^^ 

^ , 6 «QDj.m a Doin a (m -3) 


/ / (- 

o h 


0" 


2 2 

4Drem g 


-) rdrd0 


k * Cft ( -6i^D^m^+D^(tn4-3m^ ) -4D^Qin^) ( a-b) | 




r9“ 


2it a 


^21 “ / / (Oq (m^-«l)rdrd© 


0 b ^ 


CSt(Dg(i/-m2) (S^) 


= k 


12 


1^21 ■* ^12 


27t a 


‘22 


/ / D© «I (m^-1)"^ rdrd© 

Ob ? 


k 22 =* CJt Dg(ra^~l)^ ln(a/b) I 


2% a 


“2 o / 2 , \ / _2 


k 23 « / / (~2i^D^ ^ (m^-D+Dft f- (m^-l)(m ^2))rdrd0 


o b 


©^r r 


© r 


k23 = Cii(-2'^QD^(m^-l)+D£^(in^-»l)(m^~2))(a-b) j 




©' 


2-r a 

k24 = / / (-6i^^a^(m^-l)+DQCc^(m^-l)(i» -3))^^^*^® 

o b 


2 ^2 


k24 = Ck ( -6'^ gBy { ) 4-DQ(m^~-4m^ ! 
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Hi " J ) TdT < i9 

o b ^ ^ ^ 


« Cii(-2iQD^m^+DQ(ro'^-2m^)-2D^ni2) ln(a/b) 




Stc ^2 2 

^32 * ^ (ia^'*l)+ ^ (m^-l) (m^-2)rdrd0 


o b 


Cit(-2i-.D^(m2-l)+D^(m^-l)(m^-2)(a-b) 


* ^23 

I ^32 * ^23 

271 a 

^33 * •/* (4D^c^-2‘iQD^a2(m^-2)-2‘^QD^a^(m^-2) 

0 b 

+ ( ra^-2 ) ^+2D^m^p^ ) rdrdG 

■ 2.2 

^33 = C7i(4Dj.-4iQD^(m^-2)+D0(in^-2)^+2Dj^ni^)(^^^) 
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2-11 a 

= / / (12D^a^r-6'iQDj.a^r(m^-.2)-2"^ gDya^r(m^-3) 

d b 


•34 


+ D^a^rCm^-S) (m^-2) 44D^n?p^r) rdrdQ 


QDj,(m^-3)+DQ{m^-3) (10^-2) 


o 3.3 
2N/a -b 


^4D^m-)(^) 


k.^4 = Cg(l2Dy~"^QDy(8m^->18)-hDQ(m'^~5m^-h6)-f4D^m^) ; 


2^/a^-b^ 


.® ^igD^m^a^ DgmV(in^-3) 


‘41 


/ / (- 

0 b 


)rdrd6 


» Di ( -6i gD^m^+Dg( ) ( a-b) 


« k 


14 


^41 * ^14 


2it a 


k 


42 


/ / (-6'^gD^a^(m^-l)+DQa^(m^-l)(m^-3))rdrd© 


o b 


2 2 

* Crc(^6igD3^{in^-l)+Dg(!n'^-4n?+3)(-2^-5^) 


^24 


^42 * ^4 
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2tz a 

^43* / / (12D^a^x-6lQD^a^r(m^-2)-2'^QD^a^r(m^-3) 

o b 

+ D^a^r(ni^-3)(m^-2)+4D^Qin^p^r)rdrd9 

3 3 

* Cjc ( 12D^«»*^0D^( Bmp— 18 ) ■^D0( ni^— Sm^ +6) +4Dj,^n^ ) (~ — - g^* *) 


I ^43 ” ^34 I 

2it a 

^44 ^ f f (36D^a^r^-6'^0D^a^r^(in^-3)-6'^QD^a^r^(m^-3) 
o b 

+ Dgoc^r^ ( m^-3 )^+8D^n!^p^r^ ) rdrd9 

I " — — 'I . .. — P H .! ■■ .. — ,l>.. — .._ — - ' ^ 1. . I I 

. ^ 4 /j. I 

k44 =» (36Dy~12i0D^(ro ~3)-HDQ(m >3r+8D^m^)(^-^) | 

3.3 MASS MATRIX FOR ORTHOTROPIC PLATES 

Since the mass matrix is independent of material properties, 
Dq,D^ and Dj^» the expression for mass matrix of polar ortho- 
tropic plate will be same as that for isotropic plate as given 
in Section 2,2.2. 

Determination of natural frequency similarly follows as 
given in Section 2.2.3 and the results and discussion are given 
in Chapter 4. 
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CH/i?TER 4 

RESULTS AND DISCUSSION 


4.1 ISOTROPIC PLATES 

4.1.1 Isotropic Annular Plates 

The analysis presented in Section 2.2 has been applied for 
all nine possible combinations of clamped. Free and simply 
supported annular plates. Nondimens ional frequency parameters 
are obtained for all the cases and for different ratio of 
radii (b/a). For the three cases of free inner edge, the values 
are tabulated (Tables 4.1, 4.2 and 4.3) and for all the cases 
graphs are given in Figs. 4.1 to 4, g . The accuracy and con- 
vergence of the method is checked by determining natural fre- 
quencies with one, two, four and eight elements. The results 
are compared with the exact values [3 ]• From the graphs it 
can be seen that as (b/a) approaches zero each of the curves 
tends to a finite limit in all but two cases. The two exceptions 
are for the mode characterised by one nodal diameter for the 
cases of a plate with a free outside edge and either clamped or 
simply supported inner edge. The limiting frequency for each 
of these modes is zero since the motion is simply a rigid-body 
rotation of the plate about a diameter. It can also be seen 
that a variety of relationship between natural frequency and 
ratio of radii are possible. These range from the situation 



Comparison of estimates of 
with exact values; S 


NE 0.1 0.2 


1 4,91 4.78 

2 4.87 4,73 

4 4,86 4.72 

8 4.86 4.72 

Exact 4.86 - 

1 14.66 14.07 

2 14.24 14,13 

01 4 13.97 13.90 

8 13.91 13.71 

Exact 13.88 

1 29,30 

2 26.75 

0 2 4 25.67 

8 25.49 

Exact 25.45 

1 33.43 

2 24.82 

1 0 4 29.52 

8 29 .46 

Exact 29.41 

1 60.28 65.19 

2 67.60 60*62 

11 4 51,36 49.80 

8 49.07 47.19 

Exact 48.08 - 

1 97.80 78.08 

2 81.05 77.94 

4 73.76. 70.61 

8 70.23 68.23 

Exact 69.23 


27.42 

25.76 

25,06 

24.81 


34 .65 
31.61 
31.44 
31.41 


n m 


0 0 


1 


2 
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e 4.1 

frequency parameter 

■F isotropic annularplate 1-1= 0.3 


0.3 

0.4 

0.5 

4.71 

4.79 

5.09 

4.67 

4.77 

5.08 

4.67 

4.76 

5.08 

4.66 

4.77 

5.07 

4.66 

- 

5.07 

13.82 

13,66 

13.51 

13.87 

13,27 

12.66 

13.40 

12.52 

11.96 

12 .99 

12.15 

11.71 

12.82 

— 

11.62 

25.45 

23.88 

23.01 

24.77 

23.79 

22.99 

24.33 

23.43 

22.65 

24.00 

23.01 

22.29 

24.12 

— 

22.31 

41.08 

53.20 

74.43 

37.20 

47.66 

66.10 

37.06 

47.49 

65.77 

37.05 

47.47 

65 .45 

37.01 

— 

65 .76 

71.51 

79.41 

94 .30 

58.49 

64.15 

81.58 

49.31 

55.58 

72.17 

47.00 

53,84 

70.07 

45.77 

- 

69.89 

75.05 

80,92 

97.08 

73.71 

76,25 

90,76 

67,36 

70.20 

83.26 

66.41 

68.04 

80.90 

65.10 


81.13 
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Table 4.2 

Comparison of estimates of frequency parameters [wa^Cfh/D)'^*^]’ 
with exact values. C~F isotropic annular plate *= 0.3 


n 

m 

NE 

0.1 

0.2 

0.3 

0.4 

0.5 



1 

10.27 

10.52 

11.50 

13.65 

17.76 



2 

10.21 

10.44 

11.44 

13.61 

17.72 

0 

0 

4 

10.17 

10.41 

11.43 

13.60 

17 .72 



8 

10.16 

10.41 

11.42 

13.60 

17.72 



Exact 

10.16 

- 

11.42 

- 

17.68 



1 

23.66 

21.68 

21.07 

21.27 

23.06 

- 


2 

21.97 

21.63 

21.21 

21.06 

22.79 

0 

1 

4 

21.45 

21,28 

20.53 

20.26 

22.31 



3 

21.30 

20.89 

19.87 

19.82 

22.14 



Exact 

21.15 

“ 

19.50 

- 

21.98 



1 

62.52 

44.41 

36 .77 

33.30 

33.04 



2 

38 .26 

36.05 

34,10 

32.69 

32.94 

0 

2 

4 

35.31 

34,31 

33.16 

32.08 

32.49 



8 

34.77 

33.90 

32.65 

32.09 

31.67 

i' 


Exact 

34,53 

— 

32.55 

— 

32.05 



1 

51.68 

54.01 

66.99 

90.51 131.27 



2 

39.88 

43.21 

51,95 

67.46 

94.20 

1 

0 

4 

39.65 

43,09 

51.80 

67.27 

93.67 



8 

39.53 

43.03 

51.76 

67.11 

92.86 



Exact 

39.49 

— 

51,72 

“ 

93.85 



1 

79.80 

78.55 

93.15 

127.94 201.42 



2 

91.21 

92.02 

89.18 

102.72 138.51 

ii : 

1 

4 

66.74 

64.14 

65.43 

77.37 JD3.97 



8 

62.06 

59.54 

61.37 

73.53 

98.84 

III ; 

p; 

11 ' , ' 


Exact 

59-98 

- 

59.81 

— 

97.32 



1 

351.56 

304.69 

348.44 

467.26 7JE2, 96 



2 

176.58 

118.77 

104.79 

112.79 144.00 

1 

2 

4 

92,90 

87.62 

84.51 

91.18 113.88 



8 

85.58 

82.77 

81.12 

87.71 107.85 



Exact 

83.44 

- 

78,98 


107.56 
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Table 4.3 

Comparison of estimates of frequency parameters 

with exact values; F~F isotropic annular plate; 0.3 


n m 


0 2 


1 0 


1 1 


2 0 


2 1 


NE 

0.1 

0.2 

1 

5.31 

5.15 

2 

5.30 

5.15 

4 

5.31 

5.16 

8 

5.29 

4.95 

Exact 

5.30 


1 

8.86 

8.53 

2 

3.83 

8,48 

4 

8.79 

8.45 

8 

8, '78 

8.44 

Exact 

8.77 


1 

21.43 

20.51 

2 

20.67 

20.54 

4 

20.59 

20.42 

8 

20.50 

20.04 

Exact 

20,49 

- 

1 

44.39 

46.95 

2 

38 ,49 

41.82 

4 

38 .37 

41.78 

8 

38.25 

41.71 

Exact 

38.17 

- 

1 

92.14 

102.23 

2 

93.24 

77.25 

4 

62.96 

61.16 

8 

60.52 

58.18 

Exact 

58.99 

— 


0.3 

0.4 

0.5 

4.91 

4.61 

4.27 

4,91 

4.61 

4.27 

4.90 

4.62 

4.17 

4.81 

2.73 

4.00 

4.91 

- 

4.08 

8,40 

8,64 

9 .32 

8.37 

8.62 

9.32 

8.36 

8.61 

9.32 

8.36 

8.61 

9.30 

8.36 

— 

9.31 

20.32 

20.31 

20.49 

20.13 

19 .32 

18 .86 

19.41 

13 .14 

17.76 

18.76 

17.61 

17.15 

18.34 


17.18 

57.36 

76 ,04 

108.29 

50.47 

65.84 

92.46 

50,41 

65 .77 

91.96 

50.35 

65 .69 

91.71 

50.30 


91.36 

119.11 

133.17 

155 .45 

74.92 

85.45 

113.22 

62.53 

73.83 

99.17 

59.99 

71.86 

97-15 

58.82 


96.33 
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where (a) frequency constantly increases, (b) frequency initially 
decreases and then increases, and (c) frequency continuously 
decreases as ratio (b/a) increases. 

In all the cases it can be seen that the mode of vibration 
^cr-io'cl-e. 4v'<\Tm^2.l:ev 

with one nodal d-lamets r and zero nodal c-ire-jie is the fourth 

, . , , ^ , OiVc U. 

mode of vibration and that with zero nodal di r aiaetej and three 

nodal is of higher order. Hence, in Table 4.4 the 

lowest four mode of vibrations are tabulated for all the cases 

v/ith ratio of radii (b/a) = 0.1 and are compared with exact 

values. It can be seen from the table that for higher mode of 

vibration the accuracy decreases. It is because of the fact 

that in the analysis the effect of rotary inertia and additional 

deflection caused by shear forces are not taken into account. 

4.1,2 Isotropic Circular Plates 

V^hen the ratio of radii (b/a) is reduced to 0.001, the 
presence of hole is seen to have very sraall effect on the 
natural frequencies. Hence, to determine the natural frequencies 
of circular plates, annular element is used and the ratio of 
radii is taken to be 0.001. Natural frequencies are obtained 
for all the three cases of clamped, free and simply supported 
outer outer edge conditions. The results are compared with 
exact values [3 ] in Table 4.5. In order to compare the effi- 
ciency of using annular elements, the values obtained for a 
free -plate is conpared with that obtained by using a 3x12 grid 
sector element with 55 degrees of freedom [ 7 J vdiich is also 



Table 4.4 

Comparison of estimates of frequency parameters[wa^-($h/D)^*^} 
with exact values. Isotropic annular plate 0.3,(a/b)=0.1 
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Table 4.5 

Coinpar'ison of es'bima'tes of frequency paranie'ter[wa^(5h/D)^*^3 
with exact values. Isotropic circular plate : = 0,3 


Simply supported plate Clamped plate Free plate 


m 

FEM 

Exact 

FEM 

Exact 

FEM 

Exact 

Ref 

[7 3 

0 

4.94 

4.97 

10.22 

10.24 

- 

— 


1 

13.92 

13.91 

21.31 

21,25 

- 

- 

- 

2 

25.71 

25.70 

35.11 

34.80 

5.45 

5.24 

5.90 

0 

29 .70 

29.70 

39.77 

39.80 

9.00 

9.06 

8.98 

1 

49.26 

48.60 

62 .49 

60.80 

20.51 

20.50 

20.24 

2 

71.44 

70.10 

87.39 

84.60 

35.38 

35.50 

36.01 

0 

74.21 

74.10 

89.15 

89.10 

38.45 

38.40 

38.12 
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given in Table 4,5, it can be s^en that the values obtained 
by using annular element is very 

than the values obtained by using the sector element. 

4.2 ORTHOTROPIC PLATES 

4.2.1 Polar Orthotropic Annular Plates 

Natural 4) . ,^ 4 -^iined for all the nine cornbina- 

iMstural frequencies are obtar^^'^^- 

tions Of supported edge conditions. 

Frequencies calculated for 

for different hole sizes. In 4.5 to 4.8, the funda»ental 

freeuencv nt hole sizes and for different 

rrequency parameters for different 

rigidity ratios are given for S-F- 

Table 4.9 giv„ the nondimensionul f^^uency parameters for all 

the nine combinations of edge ooddltio- for b/a = 0.1. Numeri- 
cal values a« K, • f [it] «''f the three cases of 

iUBs are available in ref 

free inner edge and for rigidity . 

frequencies for those cases are ^ood agree.en is 

found. [It i, „„ted that »e values in ref. hi 

approximate values obtained by Bayleigh-Ritt method]. As a 

Chech for then t , 's the values obtained for rigidity 

tne present analysis r , t. 

ratio 1 0 i ^ t values of isotropic plates. 

' IS compared with exact _ , ^ * 0 for 

Table 4 in 4 .U +abulated for n = 1 and m 

•^0 the frequencies are t . 

of edge conditions, 

® cas03 combine tiu -hration, 

seen that ta « at higher mode of vibration, 

the accuracy decrease® 


id this is because of the i gnu 


ance of additional deflection 
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Table 4*6 

Variation of fundamental frequency parameter[u)a^(3h/p^)^‘^] 
with rigidity ratio (D^/D^) and ratio of radii (b/a)j ortho- 
tropic S”“F plate J ^xG ” 0»3b» ~ 0*3 


Rigidity 


Ratio 

of radii (b/a) 


ratio 

( Dq / Dj ,) 

0.1 

0.2 

0.3 

0.4 

0.5 

0,2 

2.29 

1.85 

1.88 

2.34 

2.77 

0.4 

3.35 

2.89 

3.03 

2 .82 

5.57 

0.6 

3.93 

3.68 

3.97 

3.61 

- 

0.8 

4.49 

4.16 

4.22 

4.45 

5.10 

1.0 

4.89 

4.89 

5.00 

4.95 

5.84 

(4.86) 


(4.66) 


(5.07) 

1.2 

5.23 

5.21 

5.04 

5.50 

6.78 

1.4 

5.49 

5.48 

5.30 

5.60 

3.04 

1.6 

5.69 

5.76 

6.08 

5.75 

7.26 

1.8 

6.03 

6.00 

6.28 

5.81 

7.65 

2.0 

6.21 

6.33 

6.42 

6.53 

5.08^ 


6.11* 




7.14'' 


Values in parentheses are exact values of isotropic plate 
^Values from Ref. [i?] 
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Table 4.7 

Variation of fundamental frequency parameter[t»ia^{^h/D^)*^*^] 
with rigidity ratio (Dq/D^) and ratio of radii (b/a)jortho- 


tropic C* 

-F plate I 

“re = ° 

.35 ; 

= 0.3 


Rigidity 
ratio , 

(Dg/Dr) 


Ratio 

of radii (b/a) 


0.1 

0.2 

0.3 

0.4 

0.5 

0.2 

7.58 

7.99 

9.68 

7.26 

10.02 

0.4 

8.61 

8.87 

10.18 

12,51 

- 

0.6 

9.35 

9.52 

11.24 

12.66 

16.54 

o.a 

9 .85 

10.03 

11.54 

13.27 

16.93 

1.0 

10.33 
( 10.16) 

10.53 

11.56 

(11.42) 

13.37 

17.91 

(17,68) 

1.2 

10.74 

10.99 

11.91 

13.88 

19.42 

1.4 

11.08 

11.33 

12.37 

14.28 

19.72 

1.6 

11.34 

11.80 

12.77 

15.33 

20.80 

1.8 

11.66 

12.04 

13.14 

14.97 

19.76 

2.0 

11.92^ 

11.45 

12.21 

13.45 

14,94 

19.86 

18.62* 


Values in parentheses are exact values of ositropic plate 
^Values from Ref. Cl7] 
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Table 4.8 

Variation of frequency parameter [wa^(Jh/D^)°*^] with rigidity 
ratio of radii (b/a). Orthotropic F-F plate; “ 0.35 ; 

= 0.3, n - 1, m » 0 


Rigidity 

ratio 

(Dq/D^) 


Ratio 

of radii (b/a) 



0.1 

0.2 

0.3 

0.4 

0-5 


0.2 

3.93 

3.27 

3.06 

3.07 

3,36 


0.4 

5.90 

5.24 

5.03 

5,10 

5.44 


0.6 

7.16 

6.58 

6.38 

6-51 

7.04 


0.8 

8.11 

7,64 

7.46 

7.66 

8.28 


1.0 

8.90 

(8.77) 

8,52 

8,39 

(8.36) 

8.64 

9.42 

(9.31) 


1.2 

9.57 

9.29 

9.23 

9.52 

10.29 


1.4 

10.17 

9.98 

9.97 

10.31 

11.12 


1.6 

10.71 

10.60 

10.65 

11.06 

12.00 


1.8 

11.21 

11.18 

11.28 

11.73 



2*0 

11 •68„ 
11.34 

11.71 

11.87 

12,39 

13.48 

13.09* 



Note * Values in parentheses are exact values of isotropic 
circular plate 


Values from Ref- Cl7 3 



Variation of frequency parameter [wa^Cf with rigidity ratio (Dq/D^) 
Orthotropic plate? - 0.35, « 0,3, (b/a) = 0.1, n = 0, m = O 
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Table 4,10 . 

Variation of frequency parameter with rigidity ratio: Orthotropic plate = 0.35, 

"^0 = 0,3 (b/a) = 0 . 1 , n = 1 , m = 0 


50 


1 ^ 

I (X. • 











CO 

\0 

r- 


l> nO 

r^ 

o 

iH 

CO 


I o 

I ^ 

%' 

NO 

¥ 

Q 

« 

tn CM 

’ • ¥ 

CD 

• 

CO 

¥ 

vO 

-¥ 

O 

¥ 

r-i 

• 

1 


\0 

NO 


r- c^ 

c- 

CD 

CO 

00 

On 

1 

CM 

CM 

CM 

CM 

CM CM 

CM 

CM 

CM 

CM 

CM 

1 

1 










iv 

1 A 











42 

O 

CO 

CM 

CM 

lO 

ON in 

16 

CO 

Nt 

CD 

00 

O 

iH 

lO 

o 

• 

» 

% 

¥ 

* • h 

« 

.¥ 

« 

¥ 

¥ 

i 


CO 



Nt O 

lO 

in 

lO 

nO 

NO 

i 

1 ^ 

0N 

On 

ON 

CTn 

ON 

CJN 

ON 

ON 

0^ 


. 










1 

lO 



CO ^ 

CO 

CO 

CD 

o 

o 

1 J 

1 o 

CO 

r- 

o 

r-i 


O 

CM 

t- 

o 

CO ■ 

i ^ 

• 


# 

• • 

• 


¥ 

• 

¥ 

1 


ON 

ON 

o 

oo 

o 

r-i 

r-i 

r-i 

CM 

1 

1 nO 

sO 

NO 

I> 

l> NO 

r- 

r- 


r- 

l>. 


j 










1 

j 










1 

ON 

I> 

o 

o 

CM 

lO 

r-1 

ON 

CO 


I CO 

' CO 

uf) 


ON 

ON 00 

nO 


r-i 

NO 

ri 

t * 

« 


• 

¥ 

* • ¥ 

* 

¥ 

m 

¥ 

¥ 

P-, 


CO 

ON 

On 

oo 

r-i 

CM 

CO 

CO 



fp 

H 

tH 

iH 

CM CM 

CM 

CM 

CM 

CM 

CM 












1 

I CM 

CO 

CM 

in 

H O 

O 

r- 

in 

in 


(0 

vO 

00 


CM 

(JnnO 


o 

nO 

-H 

O 


n 

• - 


¥ 

■ ♦ » 

¥ 

¥ 

« 

¥ 

¥ 

1 ^ ■ 

vO 

r- 

CD 

o 

ON in 

ri 

r-i 

CM 

CO 

CO 

1 

vO 

vO 

NO 

vO 

vDnO 

l> 


r- 


c^ 












1 





r--v 






1 

CO 

CO 

CD 

0N 

CO CM 

CO 

CM 

lO 


o 

1 ^ 

CM 

O 

CO 

nO 

^ r- 

CM 


i> 


r-i 

I 


* 

¥ 

¥ 

' » '• 

« 

• 

¥ 

¥ 

¥ 

CO 

r-i 

CM 

CM 

CO 

^ r-i 

in 

in 

nD 

o 

CO 

in 

lO 

lO 

in 

in tn 

in 

tO 

lO 

in 

in 

















r->. 







CD 

O 

o 

nO 

00 N;f 

r-i 

N0 

nD 

04 


tu 

CO 

0N 

CO 

Nt 



CM 

o 

CD 

nD 

1 i 

'« 


¥ 

• 

■ ♦ • 

¥ 

¥' 

¥ 

¥ 

¥ 

ix 

in 

l> 

c^ 

o 

•hcM 

CM 

CO 


xt 

tO 





r-i 

r-i r-i 

'r-i 

rH 

r4 

r-l 

r-i 






x«.r 


















tn 

co 

o 


CM in 

ON 

CM 

ON 

CO 

O 

CX 

\0 

CM 

o 

o 

0\r-i 

00 

00 

f- 

0- 

nO 

I 1 1 

% 

» 

¥ 

¥ 

' « ¥ 

• 

¥ 

¥ 

¥ 

¥ 

o 

0S 

ON 

o 

CD 

00 rHi 

00 

00 

00 

CD 

CO 




r-| 

r-i CM 

r-i 

r-| 

r-i 

r-i 

tH 












1 











1 

CM 

rH 

iD 

o 

ON CO 

r- 

tn 

o 

NO 

CNi 

^ PU 1 


lO 

nO 

00 

OCD 

r-i 

CO 

in 

nD 

CO 

■ 1 ‘1 

% 

# 

« 

• 

• • 

• 

¥ 

¥ 

¥ 

¥ 

CO 1 
I 1 

CM 

CM 

CM 

CM 

CM CO 

CO 

CO 

CO 

CO 

CO 


iH 

rH 

,—1 

r-i r-i 

r-i 

H 

»-i 

f-4 

r-l 

1 ^ 1 
1 ^ 1 











CM 


NO 

00 

o 

CM 


N0 

€0 

o 

I ® ' 1 

ii 


¥ 

¥ 

* 

• 

¥ 

¥ 

¥ 

¥ 

1 Q 1 

o 

o 

o 

o 

r-l 

•r-i- 

H 

>-i 


CM 


I 


I 

I 


o 

f—f 

a 


o 


•H 

a 

o 


u 

-p 

o 

if) 


iH 

O 


(/) 

Q) 

r~\ 

0 

> 


CO 

!l 

G 

O 

II 

c 

N 

O 

D 

> 


-p 

u 

<15 

X 

o 

<!> 

f-J 

fO 

a> 

Q 

x: 

-p 

c 

Q) 

Pi 

fO 

a 

c 

•H 

to 

m 

:3 

rH 

m 

> 



51 


caused by shear forces in the analysis. It can be seen that 
the frequency depends on rigidity ratio and hole sizes. The 
graphs given (Figs. 4.10 - 4.27) shows the variation of frequen- 
cy with rigidity ratio and variation of frequency with hole 
size. In all the cases it can be seen that the frequency in- 
creases with increase in rigidity ratio and the variation is 
linear for D^/D^ >1. In some cases, where the variation of 
frequency with rigidity ratio is small, even in the range 
Dq/D^ < 1, the variation is seen to be linear. From Fig. 4. 17 
it can be seen that the rate of increase in frequency with 
rigidity ratio is more in the range Dq/D^ < 1 and the rate of 
increase in frequency decreases with increase in rigidity 
ratio. Another thing to be noted is the variation of frequency 
with ratio of radii. In the case of a F-F plate as shown in 
Fig. 4.15, the frequency initially decreases and then increases 
with hole size, and the decreasing -tendency is found to be 
more in the range Dq/D^ < 1. For a value of rigidity ratio 
2.0, there is no decrease in frequency with increase in hole 
size, but frequency increases continuously with incacease in 
hole size. For all the cases the effect of hole size on fre- 
quency depends upon the edge conditions. 

4.2.2 Polar Orthotropic Circular Plate 


Natural frequencies are obtained for all the three cases 
of clamped, free and simply supported edge conditions. In 
order to study the variation of frequency with rigidity ratio. 




. ^ 


Frequency parameter 



Fig* 4*10 Nondimensional frequency parameter 
Vs rigidity ratio for S-F polar 
orthotropic platej o-i 
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frequencies are calculated for rigidity ratio ranging from 
0»2 to 2,0 in steps of 0,2, Since numerical values are not 
available for comparison, the accuracy is checked by conparing 
the value obtained for rigidity ratio 1 (which corresponds to 
isotropic plate) with exact value of isotropic plate. The 
fundamental frequency parameter is given in Table 4.11 and llie 
variation of fundamental frequency parameter with rigidity 
ratio is shown in Fig, 4.28, As expected the natural frequency 
increases with rigidity ratio and the variation is somewhat 
linear also. Table 4,12 gives the four lowest frequency para- 
meters of orthotropic circular plate for a given rigidity 
ratio 0,8, Here also no numerical values are available for 
comparison, but it can be seen that the values obtained are 
less than and close to the corresponding frequencies of isotro- 
pic plate as expected because the rigidity ratio is 0.8, 
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Table 4.12 
2 0 5 

Frequency parameter [wa (^h/op] * for orthotropic circular 
plates. (Dq/dP= 0,8> = 0.35 and = 0,3 


Outer edge n m Frequency 
condition parameter 
of the plate 


Free 


Clamped 


0 2 4,86 

1 0 8,44 

1 1 19.45 

1 2 32.90 

0 0 10,11 

0 1 18,99 

0 2 23.33 

1 1 58,24 

0 0 4.63 

0 1 13,48 

0 2 20.95 

1 1 44.06 


Simply supported 
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CH/:iPTER 5 

CONCLUSIONS 


5.1 CONCLUSIONS 

Vibration response of isotropic as well as polar orthotro- 
pic plates (circular and annular) are studied using finite ele- 
ment method. Annular finite elements having only four degrees 
of freedom, these being the displacements and slopes at inner 
and outer radii, are used, for the analysis. In Chapter 2, the 
vibration response of isotropic plates are studied and the 
estimates are compared with exact values. Then its applicabi- 
lity is extended to polar orthotropic plate also. In the case 
of isotropic plate, when the hole size is reduced to O.CKDl of 
outside radius, the results obtained are seen to be very close 
to the exact values of solid circular plates. Hence, the vibra- 
tion response of polar orthotropic solid circular plate is 
studied using the annular element itself. The material singula- 
rity at the centre of the circular plate, which is a major 
problem in the analysis by other methods, does not arise in the 
present analysis since here the solid plate is replaced by an 
annular plate having ratio of radii 0,001. 

In the case of orthotropic plates, it is seen that the 
frequencies depend on rigidity ratio as well as hole size. The 
frequency increases with the increase in rigidity ratio, but 
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the effect of hole size is different for different edge condi- 
tions. The efficiency of using annular element is checked by 
comparing the results with the results obtained by using a 
sector element [7 j. It is seen that the annular element with 
4 degrees of freedom gives more accurate results then the re- 
sults obtained by using sector elements having even 55 degrees 
of freedom. The results would have been better if the effect 
of rotary inertia and transverse shear deformations were 
included in the analysis. 

In fairness it must be pointed out that the use of annular 
element described here is restricted to complete annular and 
circular plate, unlike the sector element [7], which can be 
used even for segments of annular and sector of circular plates* 

5.2 SUGGESTIONS FOR FUTURE V/ORK 

(i) In the present analysis the effect of rotary inertia 
and effect of additional deflection caused by shear force were 
not considered, it is worth to study the problem including these 
effects. 

(ii) The applicability of this annular finite element can 
be extended to study vibration response of rotating polar ortho- 
tropic plates with any thickness variation in radial direction* 
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( iii) The effect of damping which is neglected in the present 
investigation tan be included in the free and forced vibration 
response analysis of polar orthotropic plates. 

(iv) An experimental study can be done to verify the results 
obtained by using the annular finite element. 
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